Doppler cooling with coherent trains of laser pulses and tunable "velocity comb 
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We explore the possibility of decelerating and Doppler cooling an ensemble of two-level atoms 
by a coherent train of short, non-overlapping laser pulses. We derive analytical expressions for 
mechanical force exerted by the train. In frequency space the force pattern reflects the underlying 
frequency comb structure. The pattern depends strongly on the ratio of the atomic lifetime to the 
repetition time between the pulses and pulse area. For example, in the limit of short lifetimes, the 
frequency-space peaks of the optical force wash out. We propose to tune the carrier-envelope offset 
frequency to follow the Doppler-shifted detuning as atoms decelerate; this leads to compression of 
atomic velocity distribution about comb teeth and results in a "velocity comb", a series of narrow 
equidistant peaks in the velocity space. 

PACS numbers: 37.10.De, 37.10.Gh, 42.50.Wk 



Laser cooling is one of the key techniques of modern 
atomic physics Radiative force originates from 

momentum transfer to atoms from a laser field and 
subsequent spontaneous emission in random directions. 
Doppler effect makes the force velocity-dependent. 




Mode-locked laser 



FIG. 1: (Color online) Schematic of a typical experimental 
setup. An atomic beam is slowed and cooled by a train of 
laser pulses. Phase and shape of pulses may be varied in time 
to attain optimal cooling. 

Here we develop a systematic theory of Doppler cool- 
ing by a coherent train of short laser pulses (see Fig.[T|). 
A qualitatively new effect comes into play: atomic 
quantum-mechanical amplitudes induced by subsequent 
pulses interfere resulting in a periodically varying radia- 
tive force as a function of frequency. This structure of 
the force reflects the comb-like pattern of Fourier image 
of the pulse train, the so-called frequency comb (EC) [J] . 
Here we derive the force and show that for sufficiently 
weak pulses and long atomic lifetimes, each tooth acts as 
if it were an independent CW laser. In the opposite limit 
of short lifetimes (short compared to the repetition time 
between pulses), we recover the force due to an isolated 
laser pulse. Earlier works on mechanical effects of FCs in- 
clude proposal involving two-photon transitions [H]. Fol- 
lowing proposal , pulse trains from mode-locked lasers 
were also used in cooling experiments [tI, Q • To the best 
of our knowledge no analytical analysis of the FCs radia- 
tive force has been attempted so far and it is presented 
here. 

Notice that over the past few years the power and spec- 



tral coverage of FCs have grown considerably. A fiber- 
laser-basjd EC with 10 W average power was demon- 
strated 9*1 and the authors argue that the technology is 
scalable above 10 kW average power. The spectral cover- 
age was expanded from optical frequencies to ultraviolet 
and to IR spectral regions [l^l. These advances pave 
the road for new applications of FCs, such as the laser 
cooling. 

As an application, we consider mapping frequency 
comb to a "velocity comb". We demonstrate that dur- 
ing pulse-train cooling, continuous velocity distributions 
gravitate toward a series of sharp (of a typical Doppler 
width of m/s for strong lines and mm/s for weak lines 
such as intercombination transition in Sr) peaks in the ve- 
locity space. "Velocity combs" could be used for studying 
velocity-dependent (e.g., shape) resonances where tradi- 
tional beam techniques with their broad velocity distri- 
butions would fail [ll'|. Moreover, since groups of atoms 
with different velocities would arrive at the target at 
different times, the experiment may be carried out "in 
parallel" for many velocities (cf. molecular fingerprint- 
ing [l3l-) Notice that the moniker "velocity comb" was 
used in a work[l3| on optical pumping with FCs; we re- 
tain this label here as a natural visual for the resulting 
velocity distribution. 

In a typical EC setup, a train of phase-coherent pulses 
is produced by multiple refiections of a single pulse in- 
jected into an optical cavity. A short pulse is outcoupled 
every roundtrip of the wavepacket inside the cavity, de- 
termining a repetition time T between subsequent pulses. 
At a fixed spatial coordinate, the electric field of the train 
may be parameterized as 



E(i) — e Ep cos(wci — 4>m) git — mT) , 



(1) 



where e is the polarization vector, Ep is the field am- 
plitude, and 4>rn is the phase shift. The frequency Wc is 
the carrier frequency and g{t) is the shape of the pulses. 
We normalize g{t) so that maxg(i) = 1, then Ep has the 
meaning of the peak amplitude. While typically pulses 
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have identical shapes and (f)m = m (p, one may want to in- 
staU an active optical element at the output of the cavity 
as in Fig.[T]that could vary the phase and the shape of the 
pulses. Also repetition time and intensity of pulses could 
be controlled by varying reflectivity of cavity mirror. 

We focus on two-level systems as these are amendable 
to analytic treatment and much insight may be gained 
from analyzing the derived expressions. Technically, we 
solve the optical Bloch equations (OBE) for density ma- 
trix elements (excited and ground state populations are 
Pee and pgg and coherences peg and pge) 

Pee = -^Pee-\' 7;{P9e^eg{z,t) - peg^ge{z,t)), (2) 



Peg = -{--iScff)peg + ^^eg{z,t){pgg- Pee), (3) 



where ^cff = (5 + kc • v is the Doppler-shifted detuning 



^eg, 



ke — 2tt/(jJc and v is the atomic veloc- 



ity). The time- and space- dependent Rabi frequency is 
flge{z, t) = V,p X]m=o g{t+z/c—mT)e^'^'^ , with the peak 
Rabi frequency ftp = -^(e|D • e\g) expressed in terms of 
the dipole matrix element. Once the OBEs are solved, 
radiative force may be determined in terms of the coher- 
ence 



-p,.lm[peg^*eg] 



(4) 



where Pr = hke is the photon recoil momentum. 

We start by observing that as long as the duration 
of the pulse is much shorter than the repetition time, 
the atomic system behaves as if it were a subject to 
a perturbation by a series of delta-function-like pulses. 
In this limit, the only relevant parameter affecting the 
quantum-mechanical time evolution is the effective area 

oo 

of the pulse = flp J g(t)dt . As an illustration we will 



we will 



use a Gaussian-shaped pulse, g{t) 
assume that Tp <C T. 

We distinguish between pre-pulse (left) and post-pulse 
(right) elements of the density matrix, e.g., {pTg)i 
{Peg)r are the values of coherences just before and just 
after the m}^ pulse. Between the pulses the dynamics is 
determined by the spontaneous decay 



Peg{t) - (pS),exp - tS,s){t ~ mT) 

Peeit)^ipTe)r^^Phlit-^T)] . 



(5) 



We neglect the spontaneous decay during the pulse, since 
for femtosecond pulses, Tp7 ^ 1. Then 



9/2 



(6) 



y, where ax,y are the Pauli 



with am — cos (j)ax — sin <p a. 
matrices. Analogs of Eqs. ()5|6p were derived earlier [14l |. 
By stacking single-pulse (O and free-evolution JS]) prop- 
agators, one may evolve a given initial p over duration of 



the entire train. In Fig. [2] we show results of such calcu- 
lation for the excited state population (atom remains at 
rest). 
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FIG. 2; Evolution of the excited state population due to in- 
teraction with a train of laser pulses. The atom is initially 
in the ground state, and it is driven by a train of pulses sep- 
arated by T = 4ns and of pulse area 9 — vr/lO. Radiative 
lifetime is 16 ns. 

Now we focus on the evaluation of the radiative force. 
The laser field is present only during the pulse, so we deal 
with a sum over instantaneous forces. The change in the 
atomic momentum due to a single pulse is 



Pr 



(7) 



i.e., a laser pulse imparts a fractional momentum kick 
equal to the difference of populations before and after 
the pulse. Since < pee < 1, the maximum momentum 
kick per pulse is equal to the recoil momentum. 

By combining Eqs. ()5|6l7p we find the radiative force. 
Time evolution of population. Fig. [21 separates into two 
regimes: initial transient phase and the quasi steady- 
state (QSS) regime when radiative-decay-induced drop 
in the population following a given pulse is fully restored 
by the subsequent pulse. Doppler cooling requires many 
scattering cycles and we focus on the QSS (or "the co- 
herent accumulation" [l3]) regime. 

In the QSS regime, pee (t) = Pee it + nT) and pre- 
and post-pulse values (p^); ^ do not depend on the 
pulse number m; we simply denote these values as 
{Ple)i,r- Then Eq.Q becomes -Ap^/pr = {ple)r ^ 
(l — e~^"^). We find (pie) using non-perturbative propa- 
gators Eqs. (|5I6[) and arrive at the fractional momentum 
kick per pulse 



-Aps sin^(6l/2)sinh(7T/2) ■ 

Pr ~ cosh(7r/2) - cos2(6'/2) cos?] 

Here the Doppler-shifted phase f] is 



(8) 



77 = ((5 + k, • v)T - cj). (9) 
Finally, the radiative force is Strain — Ap^/T. 
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FIG. 3: (Color online) Frequency dependence of the fractional 
momentum for a pulse area 9 = tv/IO. The three curves differ 
by the values of the parameter 'yT (solid curve 7T — 0.1, 
dashed -jT = 0.447 and dot-dashed -yT = 2.5). The spikes 
reflect the underlying frequency-comb teeth structure of the 
pulse train in the frequency domain. Inset: friction coefficient 
/3 as a function of phase for the dashed curve of the main 
panel. 



Its typical value (sodium atom) is 6.4 mW/cm^. Typ- 
ical FC parameters (rp = 100 fs, T = 1 ns and average 
power 1 W) translate into the power per tooth of 0.1 mW, 
thereby 1^ can be attained by focusing the laser output 
to a spot of 1.4 mm diameter. Notice that a fiber-laser- 
based FC with 10 W average power has been demon- 
strated and the authors argue that the demonstrated 
technology is scalable to 10 kW average power. With this 
new generation of combs the cross-section of the interac- 
tion region may be increased dramatically. 

Now we turn to the dynamics of slowing down and 
cooling an entire atomic ensemble, characterized by some 
velocity distribution f{v,t) (time-dependence is caused 
by radiative force). To be specific, consider a typical 
use of radiative force for slowing down an atomic beam. 
The laser pulses would impinge on the atoms (see Fig. [T]) 
countering their motion. The radiative force ([8]) depends 
on the atomic velocity via Doppler shift. As velocity is 
varied across the ensemble, the maxima of the force occur 
at discrete values of velocities {n are integers) 



= {2Trn-TS + (t))/{kcT). 



(11) 



As a function of phase t] (or frequency or velocity) , the 
force spikes at the positions of the Doppler-shifted fre- 
quency comb teeth, rjn — 2tt x n, with n being integer 
numbers (see Fig. |3l In Fig. |3]we also investigate depen- 
dence on the values of parameter jT. Let us focus on 
one of the teeth (e.g., 77 = 0). As jT is increased the 
momentum kick grows, reaches maximum, and then de- 
clines; apparently for a given 9 there is an optimal value 
oijT. By analyzing Eq. (jS]), we find this optimal value to 
be {jT)^pi = 2cosh"^ (l/cos^ (9/2)) . For example, for 

9 = 7r/10, the optimal value is ijT)^^^ « 0.447, i.e., the 
radiative lifetime is roughly twice the repetition period. 

Eq. IH]) is non-perturbative. It remains valid even for 
strong laser pulses, as long as the pulses do not overlap. 
For 9 = 2tt the force vanishes since 27r-pulse does not 
redistribute population. For 7r-pulses, the r.h.s. of Eq. ([8]) 
reduces to a heqvieTicy- independent value tanh(7r/2). 

In the limiting case of weak pulses [9^1) and fast 
repetition rates (7T <C 1), Eq. ^ reduces to 



train|^T<l,e<l 

E 



7/2 



(7/2)' + (<5 + k, • vo - ((/. + 2^n) /Tf ■ 



By comparing Eq. (fTO| with relevant CW expressions 
(see, e.g., Ref. [3i])i we arrive at the qualitative picture 
where each tooth acts as an independent CW laser, inten- 
sity of which has been reduced by the factor of ^JnTp/T. 
This factor is roughly equal to the number of teeth fitting 
inside the overall frequency envelope of the FC. 

Does a single tooth have enough power for cooling? 
The relevant parameter is the saturation intensity, I^. 



The force peaks are separated by — Vn — Xc/T in 
the velocity space. The comb may have multiple teeth 
effectively interacting with the ensemble. 

Cooling can be characterized by introducing friction 
coefHcient /3, Ftrain(f + Aw) « Ftram{v) - l3{v)Av . If 
/? > 0, there is a compression of velocity distribution 
around v. In the limiting case of 7T ^ 1 oi 9 — tt the 
force does not depend on velocity, thereby /3 — and 
while the ensemble slows down, there is no cooling. The 
friction coefficient may be derived analytically from the 
force ([5]) . We plot the dependence of /? on 77 in Fig. [31 It 
acquires the maximum value at rj — fj, 



cosjy = ^ sec^ (61/2) x 



(12) 



^Js cos4 {9/2) + cosh^{jT/2) - cosh(7T/2) 



In the CW limit, this expression leads to detuning of 
7/2 below the atomic resonance as expected. One could 
optimize /3 by varying 9 or 7T. 

As the atoms slow down, they come in and out of res- 
onances with the FC teeth, leading to periodic variation 
in the sign of /?; no cooling results due to this variation. 
To keep FC teeth in resonance with the Doppler-shifted 
atomic transition, one may vary the offset phase by 
an optical element installed at the cavity output (see 
Fig. [T] and Eq. ([TT])). If for a given velocity group ini- 
tially centered at v„ip{t = 0), the phase detuning is kept 
at fj, there will be a compression of velocity distribution 
around Vmp(t). We may satisfy this requirement by tun- 
ing the phase according to 0(i) = (5 + kcVmp{t)) T — fj. 
As Vmp{t) becomes smaller due to radiative force, the off- 
set phase needs to be reduced. Using Eq.(I5]) we find the 
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required pulse-to-pulse decrement of the phase 



pIT sin^ (6l/2)sinh(7r/2) 
hMa cosrycos^ (61/2) - cosh(7T/2) 



(13) 



We would like to emphasize an important distinction 
between earlier works [^J^ and our approach. This is 
related to the difference between slowing and cooling. 
In there is no substantial compression of velocity 

distribution about the frequency teeth. One could eas- 
ily see from figures in 7] that the spread of velocities 
about a single tooth is approximately equal to half the 
distance between neighboring teeth. There is no veloc- 
ity compression (cooling) about individual teeth because 
no frequency tuning was done in those papers. Linguisti- 
cally "velocity comb" implies narrow teeth in the velocity 
space. Eq. (jl3p prescribes how to achieve this narrow- 
ing via phase tuning. Each tooth ultimately would have 
width in the order of Doppler width or smaller. 

When the phase offset is driven according to (fT3|) . the 
entire frequency-comb structure shifts towards lower fre- 
quencies. As the teeth sweep through the velocity space, 
atomic v(t) trajectories are "snow-plowed" by teeth, ul- 
timately leading to narrow velocity spikes collected on 
the teeth. Formally, we may separate initial velocities 
into groups Vjnp{t = 0) + {fj + 27rn) /kcT < v{t = 0) < 
v,np{t = Q) + {fj + 2TT{n + l)) /kcT,n = 0,±1... The 
width of each velocity group is equal to the distance be- 
tween neighboring teeth in velocity space, Xc/T. As a 
result of "snow-plowing" , the n*'' group will be piled up 
at Vn(t) — Vmpit) + nXc/T. The final velocity spread of 
individual velocity groups will be limited by the Doppler 
temperature, = h'y/2kB- 

The formation of velocity comb is illustrated in Fig. |4] 
where we consider cooling and slowing a ID ther- 
mal beam of ^^Sr by a pulse train. The initial ve- 
locity distribution is characterized by f{v,t = 0) = 
exp(— 2^1— ), where Vmp is the most probable ve- 
locity at t = 0. In this example we use the weak 
55^^5*0 — >■ 5s5p^Pi transition with 7 = 5.3 x 10''l/s. 
Parameters of the train are T = 50 ns and 6 = 0.0019. 
At the end of the process we end up with velocity comb 
separated by 13.8 m/s and of Doppler-limited width of 
7.6 mm/s (this is comparable to the recoil limit). About 
14% of the total number of atoms is "snow-plowed" into 
the teeth within 125 /is. 

Notice that by shining two counter-propagating pulse 
trains on the atoms, one could control w-positions of ve- 
locity teeth at will, as shifting phase of one train with 
respect to the other would change the balance of two 
counter-acting radiative forces exerted by the trains. 

We demonstrated that radiative force exerted by laser 
pulse trains has unique features and expands the toolbox 
of laser cooling techniques. For example, one may engi- 
neer velocity combs that may be used for studies of nar- 
row collision resonances and thresholds 11, 151. In some 




FIG. 4: Time-evolution of velocity distribution for a thermal 
beam subjected to a coherent train of laser pulses. Pulse- 
to-pulse phase offset of the train is varied linearly in time as 
prescribed by Eq. (|13|l . Atomic and pulse train parameters 
are: 7T = 0.0026, 6 = 0.0019. The optimal phase detuning is 
?7 = 0.001. Characteristic initial temperature of the ensemble 
is 293 K. 



cases, the frequency comb may be already apart of exper- 
imental setup, e.g., in optical atomic clocks [l6|. By using 
it for cooling would reduce the number of lasers. Also the 
setup does not require Zeeman slowers, whose fields may 
be detrimental for precision measurements 17]. 
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